Light scalar field constraints from gravitational-wave observations of compact binaries 
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Scalar-tensor theories are among the simplest extensions of general relativity. In theories with 
light scalars, deviations from Einstein's theory of gravity are determined by the scalar mass m s 
and by a Brans-Dicke-like coupling parameter cjbd. We show that gravitational- wave observations 
of nonspinning neutron star-black hole binary inspirals can be used to set lower bounds on cjbd 
and upper bounds on the combination iti s /^/ubd. We estimate via a Fisher matrix analysis that 
individual observations with signal-to-noise ratio p would yield (m s / v / o7BrJ)(/9/10) < 10~ 15 , 10 -16 
and 10 -19 eV for Advanced LIGO, ET and eLISA, respectively. A statistical combination of multiple 
observations may further improve these bounds. 
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I. INTRODUCTION 

Scalar-tensor theories, in which gravity is mediated by 
a tensor field as well as a nonminimally coupled scalar 
field, are popular and simple alternatives to Einstein's 
general relativity [lj-Q. Generic scalar-tensor theories 
are of interest in cosmology, and under certain conditions 
they can be shown to be equivalent to f(R) theories [3, 
5]; they have also been investigated in connection with 
inflation and cosmological acceleration (see e.g. d,0])- 

Light, massive scalars are a generic prediction of the 
low-energy limit of many theoretical attempts to unify 
gravity with the Standard Model of particle physics. For 
example, string theory suggests the existence of light 
scalars ( "axions" ) with masses m s that could be as small 
as the Hubble scale (~ 10 -33 eV). For a coherent scalar 
field to be cosmologically distinct from cold dark mat- 
ter, or to play a quintessence-like role, it must be very 
light: I0" 33 eV < m s < 10~ 18 eV. In this "string ax- 
iverse" scenario, cosmic microwave background observa- 
tions, galaxy surveys and even astrophysical measure- 
ments of black hole spins may offer exciting experimental 
opportunities to set constraints on the masses and cou- 
plings of the scalar fields Therefore, astrophysi- 
cal and cosmological observations may provide important 
clues on the relation between gravity and the other forces. 

In this paper, building on previous work [TT|, we ex- 
plore the possibility of constraining the mass and cou- 
pling of a scalar field via gravitational- wave observations 
of compact binary inspirals with Earth-based or space- 
based interferometers, such as Advanced LIGO /Virgo 
PHI, the Einstein Telescope [H011, eLISA/NGO [H 



or Classic LISA [18| . Until recently, calculations of grav- 
itational radiation in scalar-tensor theories focused on 
massless scalar fields (see e.g. fl9l - l25l |). When m s = 0, 
deviations from general relativity can be parametrized 
in terms of a single coupling constant, that for histori- 
cal reasons is usually chosen to be the Brans-Dicke pa- 
rameter wbd- The best constraint on this parameter to 
date (ubd > <^Cass = 40, 000) comes from Cassini mea- 
surements of the Shapiro time delay (26| . For massless 
scalar-tensor theories, gravitational- wave observations of 
mixed binaries can provide constraints comparable to 
(or marginally better than) the Cassini bound p7l - [29l ]. 
Recent studies explored some of the implications that 
light, massive scalar fields may have in the context of 
gravitational-wave phenomenology. For example, reso- 
nant superradiant effects induced by light massive scalars 
may produce "floating orbits" when small compact ob- 
jects spiral into rotating black holes, leaving a distinct 
signature in gravitational waves [30l | . 

The general action of scalar-tensor theories with a sin- 
gle scalar field can be cast in the form 
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where Cm is the matter Lagrangian and collectively 
denotes the matter fields. In this paper we consider a 
scalar-tensor theory with a constant coupling function 
ui(<p) = wbd and a mass term, as in [ll|. Expanding 
the potential M(<f>) about a (cosmologically determined) 
background value of the scalar field <fi = <f>o, under the 
assumption of asymptotic flatness, we have 
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Then the equation of motion for the scalar field has the 
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schematic form 
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We note that for m s = this theory reduces to the well- 
known Brans-Dicke theory [3lj ] . 

Throughout the paper we use geometrical units G = 
c = 1, so all quantities can be expressed in (say) seconds: 
for example, the mass of the sun Mq = 4.926 x 10~ 6 s. 
We assume a standard ACDM cosmological model with 
H a = 72 kms^Mpc" 1 , fl M = 0.3 and fl A = 0.7. Given 
the scalar field mass rh s (in eV), we find it convenient to 
define a quantity m s = rh s / ft with dimensions of inverse 
length (or inverse time, since c = 1), as this is the quan- 
tity that would appear in the fiat-space Klein-Gordon 



:)<y5 = 0. To convert between m s and 
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m s , it is sufficient to note that ft = 6.582 x 10~ 16 eVs. 
We will use rh s and m s interchangeably in the rest of the 
paper; the units should be clear from the context. 

Our study is the first to explore bounds on massive 
scalar tensor theories using gravitational-wave observa- 
tions of compact binaries. We only consider quasicir- 
cular, nonspinning binaries for two reasons: (1) resonant 
effects in the context of gravitational- wave detection were 
considered in [32[ , and they would produce large enough 
dephasing to prevent detection using general relativis- 
tic templates; (2) in the absence of resonant effects, the 
introduction of aligned spins would increase parameter 
measurement errors by about one order of magnitude. 
However, spin precession and eccentricity would reduce 
the errors by a comparable amount [29|, [33| . For these 
reasons, our bounds for quasicircular, nonspinning bina- 
ries should be comparable to bounds resulting from more 
realistic (and complex) waveform models, as long as res- 
onant effects don't play a role. 

In this paper we will focus on quasicircular neutron 
star-black hole binaries. The main reason is that, accord- 
ing to our previous investigation (Tlj . only mixed bina- 
ries (i.e., binaries whose members have different gravita- 
tional binding energies) can produce significant amounts 
of scalar gravitational radiation. Under the assumption 
of asymptotic flatness, dipole radiation is produced due 
to violations of the strong equivalence principle when the 
binary members have unequal "sensitivities": s% ^ S2- 
The sensitivities are related to the gravitational binding 
energies of each binary member (labeled by a = 1,2): 
s a = 1/2 for a Schwarzschild black hole, s a ~ 0.2 for a 
neutron star, and s a ~ 10~ 4 for a white dwarf [2(J [Hj]. 
Roughly speaking, dipole radiation is produced when the 
system's center of mass is offset with respect to the cen- 
ter of inertia, so mixed and eccentric binaries are the best 
observational targets to constrain scalar-tensor theories. 
A second reason to consider neutron star-black hole bi- 
naries is that dipolar radiation should not be emitted in 



black hole-black hole systems because of the no-hair theo- 
rem, i.e. the fact that black hole solutions in scalar-tensor 
theories are the same as in GR (see [1^] and references 
ther ein) . Recently, building on earlier work by Jacob- 
son [3g, Horbatsch and Burgess pointed out that scalar 
fields that vary on cosmological timescales may violate 
the no-hair theorem, so that even black hole-black hole 
binaries may produce dipolar radiation 37 1. The exis- 



tence of dipolar radiation in extended theories of gravity 
and the investigation of possible bounds on dipolar radi- 
ation are active research topics |38l - l4"oT | . 

In summary, here we will make the conservative as- 
sumption that only mixed binaries generate dipolar ra- 
diation, and we will investigate bounds on the mass and 
coupling of the scalar field coming from gravitational- 
wave observations of black hole-neutron star binaries. 
For simplicity we will set the mass of the neutron star 
to be Mns = 1.4 Mq, and we will focus on nonrotating 
black holes. This rules out by construction the possibility 
of floating orbits of the kind studied in [3(1 [Hj . 

We will consider the bounds that could be obtained us- 
ing different gravitational- wave detectors. For Advanced 
LIGO we use the fit to the expected power spectral den- 
sity given in Table I of jU (henceforth "AdLIGO" ) as 
well as the zero-detuning, high power configuration, as 
given analytically in Eq. (4.7) of g3 ( "AdLIGO ZDHP" ) . 
In both cases we assume the power spectral density to be 
infinite below a seismic noise cutoff frequency of 20 Hz. 
For the Einstein Telescope we use the analytical fit pre- 
sented in Table I of [41(, assuming a lower cutoff fre- 
quency of 10 Hz. For Classic LISA we use the analytical 
Barack-Cutler expression 43 1, as corrected in (28[. For 



eLISA/NGO we use the noise model (inclusive of galactic 
background noise) discussed in [TtJ ■ 

In general, the gravitational radiation from a binary 
in massive scalar-tensor theories depends on the scalar 
field mass m s and on the coupling constant £ ~ 1/cjbd 
(see [ll[ and Eq. Q below for a more precise defini- 
tion). Matched-filtering detection of gravitational radi- 
ation from compact binaries relies primarily on an ac- 
curate knowledge of the wave phasing. In any theory 
that differs from general relativity, the gravitational- wave 
phasing will be modified by a slight amount that de- 
pends on the parameters of the theory. In massless 
scalar-tensor theories the dominant (dipolar) contribu- 
tion scales like 1 /wbd ■ If we measure a waveform consis- 
tent with general relativity we can set upper bounds on 
deviations of the dipolar term from zero, and therefore 
we can set lower limits on wbd (27l - [29j . By computing 
the gravitational-wave phase in the stationary-phase ap- 
proximation (SPA; cf. Eq. (|16p below) we will see that 
the scalar mass always contributes to the phase in the 
combination ~ m^/weD, so that gravitational- wave 
observations of nonspinning, quasicircular inspirals can 
only set upper limits on m s / ^lobd- For large signal-to- 
noise ratio (SNR) p, the constraint is inversely propor- 
tional to p. It is useful to plot constraints that result 
from observations with SNR p = 10: this corresponds to 
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a minimal threshold for detection of the corresponding 
binary system, and in this sense it places an upper limit 
on the achievable bounds on m s (and a lower limit on 
the achievable bounds on wbd)- 
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FIG. 2. Bounds from eLISA/NGO at SNR p = 10, compared 
to Solar System and binary pulsar bounds. 



FIG. 1. Bounds on m s /y^irJ with AdLIGO, AdLIGO ZDHP, 
ET, Classic LISA and eLISA/NGO at fixed SNR p = 10. 

In Fig.[Tjwe plot the upper bound on m s / ^/lu-qti result- 
ing from neutron star-black hole binary observations with 
SNR p = 10, as a function of the black hole mass Mbh, 
for different detectors. The order of magnitude of these 
bounds is essentially set by the lowest frequency accessi- 
ble to each detector, and it can be understood by noting 
that the scalar mass and gravitational-wave frequency 
are related (on dimensional grounds) by m s (eV) = 6.6 x 
lCT 16 /(Hz), or equivalents /(Hz) = 1.5 x 10 15 m s (eV). 
For eLISA, the lower cutoff frequency (imposed by ac- 
celeration noise) / cu t ~ 10 -5 Hz corresponds to a scalar 
of mass m s ~ 6.6 x 10~ 21 eV. For Earth-based detectors 
the typical seismic cutoff frequency is ~ 10 Hz, corre- 
sponding to / cu t ~ 6.6 x 10~ 16 eV. These lower cutoff 
frequencies set the order of magnitude of scalar masses 
probed by space-based and Earth-based detectors. 

The best bounds are obtained from the intermediate 
mass-ratio inspiral of a neutron star into a black hole of 
mass Mbh ^ 10 3 M Q , as observed by a space-based in- 
strument such as eLISA or Classic LISA. In summary, we 
conclude that the most competitive bounds would come 
from space-based gravitational wave detectors, and that 
they would be of the order 
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In Fig. [2] we plot current bounds on wbd as a func- 
tion of m s coming from (i) Cassini measurements of the 
Shapiro time delay (solid black line), (ii) Lunar Laser 
Ranging bounds on the Nordtvedt effect (dotted blue 
line), and (iii) measurements of the orbital period deriva- 
tive of two binary pulsar systems. We refer the reader to 
fj~TI | for a detailed derivation of these bounds. 



On top of the existing bounds we plot, for comparison, 
future bounds from eLISA observations of neutron star- 
black hole binaries where the black hole has mass Mbh = 
300 M Q (solid lines) and M BH = 3 x 10 4 M Q (dashed 
lines). In both cases the bounds refer to observations 
with SNR p = 10. 

For each system, as discussed above, gravitational 
wave observations provide two constraints: a lower limit 
on cjbd (corresponding to horizontal lines in the plot) 
and an upper limit on m s / Jwbd (corresponding to the 
straight diagonal lines). For each of the two neutron 
star-black hole systems, a gravitational-wave observation 
would exclude the complement of a trapezoidal region 
on the top left of the plot. We see that bounds from 
space-based gravitational- wave observations will improve 
over the Cassini, Lunar Laser Ranging and binary pulsar 
bounds bound only for neutron star-black hole systems 
with Mbh ^ 10 3 M@. In a sense, these predictions are 
conservative: for example, for a gravitational-wave ob- 
servation with SNR p = 100 both bounds (the horizontal 
and diagonal lines) would improve by an order of magni- 
tude. Therefore, a single high-SNR observation (or the 
statistical combination of several observations, see e.g. 
44 1) may yield bounds on massive scalar-tensor theories 
that are significantly better than Solar System or binary 
pulsar bounds. 

The plan of the paper is as follows. In Section [TT] we 
compute the gravitational- wave phase for massive scalar- 
tensor theories in the stationary phase approximation. In 
Section HTI1 we use these results to compute the Fisher in- 
formation matrix in these theories. Section IIVI presents 
the bounds obtainable with Earth- and space-based in- 
terferometers. Section [V] discusses possible directions for 
future research. 
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II. STATIONARY PHASE APPROXIMATION 
IN MASSIVE SCALAR-TENSOR THEORIES 

As discussed in SectionU we focus on scalar-tensor the- 
ories with constant coupling and a massive scalar field, 
following the treatment and notation of ll| . We consider 



where 



the inspiral of a binary system composed of two compact 
objects with masses mi and mi- For consistency with 
the notation of [28( . here and below all masses are 
measured in the detector frame; they are related to the 
masses to| ' in the source frame by = (1 + z)mf\ 
Let to = toi + to 2 be the total mass of the binary, 
u = Tjm = (TO1TO2)/ (mi + TO2) the reduced mass (with i] 
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the symmetric mass ratio), and A4 

the so-called chirp masfl We will denote by / the 
gravitational-wave frequency of the binary. 

The dominant contributions to gravitational radiation 
in scalar-tensor theories are dipolar and quadrupolar. 
In the adiabatic approximation, the time derivative of 
the gravitational-wave frequency (denoted by an over- 
dot) reads [ll| 
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and 8(x) is the Heaviside function. 

The gravitational-wave and orbital frequencies are re- 
lated by / = 2/ orb , where (see e.g. Eq. (2.25c) in [20j]) 
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Following [ia |46|], we can compute the phase ip(f) of 
the gravitational waveform as a function of the wave fre- 
quency / in the stationary phase approximation (SPA): 



4>(f) = 2nft(f) - 4>(f) - 



(9) 



1 Notice that Refs. [53 . [2?| introduce a slightly different definition 
of the chirp mass in the context of scalar- tensor theories. 
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For light scalars (m s < 10 -19 eV) we know that cjbd 
must be larger than the Cassini bound wcass = 40, 000 
(see Fig. [2]), therefore we arc justified in assuming that 
{ < 1 and we can linearize in £. For scalar masses 
m s > 10~ 18 cV, strong couplings ~ 1 or larger) are 
not experimentally ruled out; however our focus here is 
on constraining light scalars of cosmological relevance. A 
straightforward calculation linearizing in £ yields 
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where dots denote contributions of higher post- 
Newtonian order. Let us define the quantity 

(12) 
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This quantity is dimensionless, since (as discussed above) 
the scalar mass m s has dimensions of inverse length (or 
time), and in geometrical units the total mass of the bi- 
nary system has dimensions of length (or time) . 

Let m2o = 10~ 20 eV be the typical scalar mass below 
which we get bounds on the scalar coupling from Solar 
System experiments [Tl[ [4?| , and let Mq be the mass of 
the Sun. We get the scaling 

, = 5.60xW^)7-^Y, (13) 
\m 2 oJ \MqJ 

so that a binary with m = 10 6 M Q (a typical target for 
eLISA) could have at most v ~ 5.60 x 10~ 9 when m s = 
in-2{). With this definition we can write 



/ 



96 

~ 

12 

5gS 2 
48 



1 



^8/3^5/3^1/3 1 * i{ai + S2 _ 2sifi2) _ ^ 



9(27t/-to s ) 



2(^to/) 2 
1 



16(ttto/) 4 



(ttto/) 2 / 3 (ttto/) 8 / 3 



6(tt/ 



(14) 



Furthermore, we note that (dubbing cither /q = to s /(27t) 
or /o = to s /7t) 

f df (f) n ®(f - /„) = e(/ - /o) f f dfur 
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where the constant can be absorbed in "time at coales- 
cence" and "phase at coalescence" integration constants 
(t c , <j) c ). The final result for the SPA phase reads 

W) = ^ " " 7 + 128( JL/)W x 
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Adding higher-order corrections in the post-Newtonian 
velocity parameter v — { / nmf) 1 /' i = (7rA / f/) 1 / 3 ?7 _1/ ' 5 to 
the standard GR phase is trivial. Eq. (2.2) of [28| lists all 
contributions up to 2PN (including spin-orbit and spin- 
spin interactions). Here we ignore spins, but we do in- 
clude all nonspinning contributions up to 3.5PN order: 
cf. Eq. (3.3) and (3.4) of @. 

Arun recently argued on general grounds that dipolar 
emission in generic extensions of general relativity should 
always introduce a term proportional to v~ 2 in the SPA 
[39l| . Our calculation shows that a nonzero mass of the 
scalar introduces additional structure in the waveform: 
in particular, Eq. (| 10[) contains terms proportional to 
v (and hence to m 2 ) that scale like v~ s and v -6 , and 
terms proportional to v 2 that scale like v~ 12 . Aran's 
argument to constrain dipolar radiation J3||, while cor- 
rect when m g = 0, is not general enough to encompass 
all scalar-tensor theories (let alone theories whose action 
contains quadratic or higher-order terms in the curva- 
ture, such as Gauss-Bonnet or Chcrn-Simons modified 
gravity [H, H(| ) . While these large negative powers of v 
could in principle be strongly dominant over "standard" 
quadrupolar radiation at small frequencies, the presence 
of the Heavisidc functions protects the SPA phase from 
a possible "infrared divergence" in the limit v 0. 

Arun pointed out that scalar-tensor theories cannot 
be easily incorporated within the original parametrized 
post-Einsteinian framework because of additional terms 
that should appear as amplitude corrections to the wave- 
form [Hj . This issue was recently addressed in jioj ]. 
where the parametrized post-Einsteinian framework was 
extended to allow for amplitude corrections. The prob- 
lem here is of a different nature, since we work in the re- 
stricted post-Newtonian approximation and therefore we 
ignore amplitude corrections "by construction". What 



happens instead is that some phase corrections are multi- 
plied by Heaviside functions 0(tt/— m s ) and 0(2Trf—m s ) 
that reflect the dipolar or quadrupolar nature of the var- 
ious contributions to the flux. This is a conceptual diffi- 
culty with power-law parametrizations of deviations from 
general relativity, but it is not an obstacle in practice, be- 
cause we know from Solar System observations that m e 
must be small. Therefore in our parameter estimation 
calculations we will assume that m s = 0, and look for 
deviations of m s from zero. Under this approximation 
the condition / > m s /7r is always satisfied, and all Heav- 
iside functions can be set equal to one. 

In massive scalar-tensor theories, corrections to gen- 
eral relativity depend on two small parameters, v and 
£. The structure of the phasing function (|16p has a def- 
inite hierarchy when seen as a multivariate Taylor se- 
ries in £ and v. However, the term proportional to v~ 12 
(which is proportional to £v 2 and formally of "third or- 
der" in the small parameters) will produce sizeable cor- 
rections whenever v ~ v 6 . More generally, it is possi- 
ble that some extensions of general relativity will con- 
tain additive corrections to the phasing of the form (say) 
av~ a + /3v~ b , with a < 1, /3 < 1, a > and b > 0. The 
parametrized post-Einsteinian framework relies on iden- 
tifying a leading-order correction to a general relativistic 
waveform. In such cases it may be hard to identify a sin- 
gle leading-order term: a simple power counting in v is 
not sufficient, unless we have independent means of con- 
straining the relative magnitude of the small parameters 
a and j3. 



III. FISHER MATRIX IN MASSIVE 
SCALAR-TENSOR THEORIES 



Following the notation and terminology of [28|, we 
work in the angle-averaged approximation and we 
write the gravitational waveform in the restricted post- 
Newtonian approximation (where only the leading-order 
term is kept in the wave amplitude) as 



Kf) = 7Ar 7/6 e W) 
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Here A = "^^573 Dl is the luminosity distance to 
the source and 7 is a geometrical correction factor: 7 = 
1 for Earth-based detectors, and 7 = V3/2 for Classic 
LISA and eLISA (cf. (28j). When the binary members 
are nonspinning the gravitational waveform depends on 
seven parameters: {ln„4, t c , <f) c , In Ai, In 77, £, v\ . 

An estimate of the accuracy in determining these pa- 
rameters can be obtained using the Fisher matrix for- 
malism. The calculation of the Fisher matrix requires 
explicit expressions for the derivative of the waveforms 
with respect to the various parameters. We refer the 
reader to [2^, |4(| for a detailed overview of the Fisher 
matrix calculation in the present context, and to Ap- 
pendix [S] for some subtleties in estimating the error on 
m s given the error on v. For theories involving a massive 
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scalar, Eqs. (2.13) and (2.14) of [28| must be modified as 
follows. Eq. (2.13f) becomes 
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survive as 

m s — > (or equivalently as v — > 0), and they correspond 
to the and K§ terms in [28j . In the absence of spins 
our waveform now depends on seven parameters (rather 
than six, as in the massless case), so we have two equa- 
tions replacing (2.13e) of [28|: 
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We remark again that we are interested in deviations 
of m s and £ from zero. Following the procedure outlined 
after Eq. (2.15) in (28|, when computing the Fisher ma- 
trix and its inverse we set m s = v = in the derivatives 
listed above (of course, all Heaviside functions are then 
equal to one). We verified the validity of this assump- 
tion by repeating the calculation with m s = 10 -20 eV (a 
mass for which we know from Solar System constraints 
that the coupling must be small [HI): our bounds are 
essentially unaffected by this modification. However it 



should be remarked that it is not possible to set £ = 0. 
The reason is apparent from Eq. (|22[) : ^ is proportional 
to £, and therefore the Fisher matrix becomes singular in 
the limit £ — > (wbd — > oo). 

By inspection we see that when v contributes to the 
SPA phase in Eq. (|T6|) it is always multiplied by the cou- 
pling parameter £, so we can only constrain the product 
£i/, or in other words the combination Sv/ujbd, where 
5v measures variations of v from zero. In practice we 
can estimate the accuracy in estimating 8v by setting 
£ = (2 + wbd) -1 equal to some small but nonzero value 
compatible with Solar System bounds. For the calcu- 
lations in this paper we chose wbd = 4 x 10 5 , but we 
verified that the measurement error on 8u scales linearly 
with l/£ for small coupling parameters. Therefore we 
present our results in terms of the combination 
which is independent of the assumed value of wbd in the 
small-coupling limit. 

Appendix A shows that once we know the statistical 
error on v, say ov, the statistical error a mg on m s can be 
estimated as 



1/2 



(23) 
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FIG. 3. SNR for observations of neutron star-black hole bi- 
naries as a function of the black hole mass Mbh for various 
detectors (AdLIGO, AdLIGO ZDHP, ET, Classic LISA and 
eLISA/NGO) at luminosity distance D L = 200 Mpc. 



IV. BOUNDS WITH EARTH- AND 
SPACE-BASED INTERFEROMETERS 

For each of the five detectors (AdLIGO, AdLIGO 
ZDHP, the Einstein Telescope, Classic LISA and eLISA) 
we present "conservative" estimates of the bounds by 
considering binaries with fixed SNR p = 10. There are 
two reasons for this choice: one is conceptual (Fisher ma- 
trix calculations are unreliable for black hole masses such 
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FIG. 4. Bound on wbd obtained and considering a six-parameter Fisher matrix at fixed SNR p — 10. Left: we set m a oc v = 
and consider a six-parameter Fisher matrix. Right: same, but for a seven-parameter Fisher matrix with f 



that p < 10) and one is practical (using p = 10 allows 
direct comparison with previous work, e.g. f2c|). 

We consider neutron star-black hole binaries where 
the neutron star has mass Mns = 1.4Mq, and we vary 
the black hole mass in a range depending on the op- 
timal sensitivity window of each detector. Bounds are 
inversely proportional to the SNR, which in turn is in- 
versely proportional to the luminosity distance of the 
binary (p ~ I/Dl). To facilitate rescaling of our re- 
sults, in Fig. [3] wc plot the SNR of neutron star-black 
hole binaries at luminosity distance Dl = 200 Mpc as 
a function of the black hole mass Mbh, for all five de- 
tectorfl Note that AdLIGO and AdLIGO ZDHP are 
very similar in terms of SNR. The same applies to Clas- 
sic LISA and eLISA: in fact, eLISA has slightly larger 
SNR when M B h < 2000 M©, mainly because the eLISA 
"armlcngth" is a factor of five smaller with respect to 
Classic LISA (see [13). 

The bounds that can be placed on m s / y/uJBD were dis- 
cussed in SectionQ] In Fig.|4]we complement those results 
by plotting the bounds on wbd , estimated as in [28| . The 
left panel shows the bounds on wbd that we would obtain 
if we considered massless scalar tensor theories, as in {28j . 
In the right panel we show that when m s ^ the bounds 
get worse by about one order of magnitude. This is ex- 
pected: we are adding one additional, highly correlated 
parameter to the waveform, and this reduces parame- 
ter estimation accuracy on all intrinsic parameters of the 
binary by roughly one order of magnitude. This degra- 
dation of the bounds is analogous to what happens when 
we add spin-orbit terms to the SPA: cf. the discussion 



2 The neutron star-black hole binaries of interest here have large 
mass ratio. For large mass-ratio systems the SNR is proportional 
to Mns/Mbh- Bounds on the scalar mass and on ojbd scale 
linearly with the SNR, so changing the mass of the neutron star 
has the trivial effect of rescaling the bounds by a constant factor 
which is very close to unity. 



around Table IV of |28|. As we could expect from the 
SNR plots, AdLIGO bounds are very similar to AdLIGO 
ZDHP bounds, and eLISA does slightly better than LISA 
for small black hole masses. Whether we consider mas- 
sive or massless scalar-tensor theories, the best bounds 
(competitive with the Cassini bound) would come from 
observations of the intermediate mass-ratio inspiral of a 
neutron star into a black hole of mass Mbh ^ 10 3 M Q , 
as observed by a space-based instrument such as eLISA 
or Classic LISA. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we have studied the bounds on mas- 
sive scalar tensor theories with constant coupling from 
gravitational-wave observations of quasicircular, non- 
spinning neutron star-black hole binary inspirals in the 
restricted post-Newtonian approximation. We found 
that neutron star-black hole systems will yield bounds 
(m s / y /UJ B ^)(p/l0) < 10~ 15 , 10'~ 16 and 10" 19 eV for Ad- 
vanced LIGO, ET and Classic LISA/cLISA, respectively 
We also found that the best bounds on wbd would come 
from space-based observations of the intermediate mass- 
ratio inspiral of a neutron star into a black hole of mass 
Mbh < 10 3 M . 

It would be interesting to drop the restricted post- 
Newtonian approximation and to consider amplitude cor- 
rections in the context of massive scalar-tensor theories 
(cf . [1^, H3] for recent work in this direction) . Further- 
more, in our analysis we have neglected time delay effects 
which arise because the massive scalar modes propagate 
slower than the massless tensor modes; while presumably 
small, it is worthwhile to investigate how these effects 
would change our bounds. Another obvious extension 
of this work would be to include spin precession, orbital 
eccentricity and the mcrgcr/ringdown waveform. All of 
these effects could improve our conservative estimate of 
the bounds on m s and wbd- Finally, it would be inter- 
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csting to estimate the improvement on the bounds that 
would result from observing several neutron star-black 
hole systems with one or more detectors (see e.g. pl.l49j 
for preliminary studies in a slightly different context). 
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Appendix A: Some subtleties in the calculation of 
parameter estimation errors 

The inversion of the Fisher matrix yields the statistical 
errors on the seven parameters chosen to parametrize the 
waveform, namely {ln.4, t c , <j) c , In M., In 77, £, i/}, and the 
associated correlation coefficients. Our goal is to get a 
bound on m s , assuming that the background value of m s 
is zero. 

In order to evaluate the statistical error of the scalar 
field mass m s , we first find the statistical error of the 
variable 



ml = 



6/5 



(Al) 



as a function of the statistical errors of v, M, ij. 

The error on a parameter 9 a can be obtained by in- 
verting the Fisher matrix 



dh 



dh 
d¥ 



(A2) 



If we consider the logarithm of a parameter, say 9 a 
ln# a , as we do here for M. and 77, then obviously 



Tab — 



dh 
din 9 



dh \ 



So switching to the logarithm of a parameter corresponds 
to a rescaling of the corresponding matrix element by 9 a . 
For the elements of the variancc-covariance matrix E ah 
we have 



-lp-1 

1 ab 



(A4) 



The transformation of the variance-covariance matrix un- 
der a change of variables is given (e.g.) in Appendix A 
of [5(|. Suppose we are given functional relations Vi{xj), 
and introduce the Jacobian D % J = 5^. Following the 

x dxj & 

notation of that paper, the result reads 

E a = (D x )- 1 S y ((D a ) T )- 1 . (A5) 
For reasons explained in Appendix A of (Hoj ]. by consid- 
ering the diagonal components of the error matrix in the 
new variables y we get a conservative estimate of the er- 
rors. In conclusion, if x% = (v , In M. , In 77) we can write 



^12/5 



dm 2 s dm 2 s £ 
dxi dxj 1 



(A6) 



M 4 



3G 



12 24 2 



Because we assume that v 
expression reduces to 



in the background, this 



(A7) 



Following common use in the gravitational-wave litera- 
ture, we will ignore subtle issues associated with the fact 
that m s > and simply use a ms = ^/er m 2 to obtain the 
conservative estimate given in Eq. (|23[) of the main text: 
if there is a 68% probability of finding m 2 s < a m z , then 



(A3) there is the same probability of finding m s < Jo^ 
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